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A microscopic expression of the residual electrical resistivity tensor is derived in linear response 
theory for Rashba electrons scattering at a magnetic impurity with cylindrical or non-cylindrical po¬ 
tential. The behavior of the longitudinal and transversal residual resistivity is obtained analytically 
and computed for an Fe impurity at the Au(lll) surface. We studied the evolution of the resistivity 
tensor elements as function of the Rashba spin-orbit strength and the magnetization direction of 
the impurity. We found that the absolute values of longitudinal resistivity reduces with increasing 
spin-orbit strength of the substrate and that the scattering of the conduction electrons at magnetic 
impurities with magnetic moments pointing in directions not perpendicular to the surface plane 
produce a planar Hall effect and an anisotropic magnetoresistance even if the impurity carries no 
spin-orbit interaction. Functional forms are provided describing the anisotropy of the planar Hall 
effect and the anisotropic magnetoresistance with respect to the direction of the impurity moment. 
In the limit of no spin-orbit interaction and a non-magnetic impurity of cylindrical symmetry, the 
expression of the residual resistivity of a two-dimensional electron-gas has the same simplicity and 
form as for the three-dimensional electron gas [J. Friedel, NuovoCimento Suppl. 7 , 287 (1958)] and 
can also be expressed in terms of scattering phase shifts. 

PACS numbers: 72.10.Bg, 72.10.Fk, 72.20.My 


I. INTRODUCTION 

The electron transport in metals and semiconductors 
is an important field of research since it crucially infiu- 
ences the efficiency, the power consumption, the size and 
the lifetime of electronic components. The use of the spin 
degree of freedom in addition to the charge is expected 
to boost microelectronics by adding new functionalities 
to existing devices. An important building block relating 
charge to spin currents is the spin-orbit (SO) interaction 
that appears also in terms of the Rashba effeclP^ in sur¬ 
face and interface states of heavy metals or in semicon¬ 
ductors in contact to those. The Rashba effect arises 
from the SO coupling in an environment with a lack 
of space inversion symmetry such as interfaces and sur¬ 
faces where Bloch-momentum-dependent spin splittings, 
known as the Rashba spin splitting, are observed in the 
band structure. The Rashba spin splitting at surfaces 
was first observed for the two-dimensional (2D) Shockley 
surface state of the Au(ll l) sn rface.l^M^ It was in vesti- 
gated for a number of clearP^ and alloyed surface^^^, 
adsorbed surface layer^HE^ and surface states of semi¬ 
conductors in contact to heavy metals.fl^ For some of 
those systems large spin-splittings had been observed. 
By this, the Rashba splitting makes possible the efficient 
application of the Edelstein effeclP^ mediated creation of 
large lateral spin-polarization. 

The scattering of Rashba electrons at impurities is 
a source of magneto-transport properties. In partic¬ 
ular we expect contributions to the planar Hall effect 


(PHE), the anomalous Hall effect (AHE) as well as the 
anisotropic magnetoresistance (AMR). Both the PHE 
and AHE ar e obse rved as a voltage transverse to the ap¬ 
plied currenlf^^^ in contrast to the AMR, which is mea¬ 
sured in the longitudinal geometry. Consequently PHE 
and AHE are characterized by the transverse resistivity 
pxy^ while the AMR is characterized by the longitudinal 
one, pxx- For PHE the in-plane component of the mag¬ 
netization M|| with respect to the interface plane and for 
AHE the out-of-plane component matter. Although 
the AMR is known since 1856^ the PHE was discov¬ 
ered more than a century later.^ Only recently it could 
be established that the PHE can originate from the AMR 
without contributions from the aheP 

Since the early nineties, with the advent of low- 
temperature scanning tunnelling microscopy (STM) the 
investigation of the scattering of surface states on im¬ 
purities proved very fruitful for studying the quantum 
behavior of matter. Crommie et alW^ reported the direct 
observation of standing-wave patterns due to the scatter¬ 
ing of the 2D electron gas off a single Ee impurity on the 
Au(lll) surface. Recently Lounis et refined this 
picture by including the Rashba effect in the descrip¬ 
tion of the 2D electron gas provided by the Shockley 
surface state of the Au(lll) surface and found a mag¬ 
netic adatom induced skyrmion-like spin texture in the 
standing wave pattern. 

In this paper we study the effect of the SO interac¬ 
tion on the residual resistivity of 2D surface or inter¬ 
face Rashba states induced by a single non-magnetic or 
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magnetic impurity and their contribution to the various 
magneto-transport properties. The impurity may not 
necessarily be a single atom. It can be any defect whose 
potential perturbation is localized. For example, it can 
be an adsorbed molecule made out of several atoms. An 
objective of this paper is to formulate a residual resistive 
tensor with longitudinal and transverse resistivity contri¬ 
butions as an extension of a recently derived microscopic 
linear response expression of the surface resistivity for a 
semi-infinite jellium modeP^^ to the 2D Rashba elec¬ 
tron gas. The impurity in our model is described by a 
scattering matrix and we consider in general that the 
potential perturbation induced by the impurity does not 
need to fulfill any symmetry such as cylindrical symme¬ 
try with respect to the impurity position or a magnetic 
moment perpendicular to the surface. Moreover, the SO 
interaction can be incorporated at the level of the im¬ 
purity besides the surrounding electron gas (see RefP). 
In that case the SO coupling contribution to the scatter¬ 
ing matrix can be added. For our application related to a 
single magnetic atom adsorbed on Au(lll), we will, how¬ 
ever, neglect this SO term when evaluating the residual 
resistivity tensor. While our focus is on the Rashba elec¬ 
tron gas, we note that recently an interesting work was 
performed on the impact of impurities on the transport 
properties of 3D topological insulators surface^^. 

We found that in absence of SO interaction and of the 
magnetism at the impurity, the diagonal contributions of 
the resistivity tensor induced by a single adatom can be 
expressed in terms of scattering phase shifts just alike the 
well-known expression by FriedeP^, 

p = Af + 1) sin^[(5^+i(5F) - , (1) 

^>0 

where the residual resistivity p, at T = 0 K induced by a 
non-magnetic impurity with a spherical potential in a 3D 
degenerate free-electron gas is related to the momentum- 
transfer cross section of electrons at the Fermi surface 
given by a sum over the orbital momentum quantum 
numbers i. Here Af is a constant prefactor proportional 
to the inverse of the wave vector k-p at the Fermi energy 
5 f of the host conduction electrons described by the free 
3D electron gas. 61 ( 6 ^) is the phase shift between the 
scattered wave function at the impurity and the unscat¬ 
tered host wave function. 

As an application we have considered an Fe impurity 
on the Au(lll) surface, both non-magnetic and magnetic 
with perpendicular and arbitrary direction of the impu¬ 
rity moment with respect to the surface plane. If the 
impurity moment has a component parallel to the sur¬ 
face plane the scattering matrix as well as the transverse 
components of the resistivity tensor are non-zero even 
when the impurity potential has cylindrical symmetry 
and a PHE is found. We follow the longitudinal resistiv¬ 
ity and the AMR as function of the spin-orbit strength 
and the transverse resistivity and the PHE as function 
of the orientation of the magnetic impurity moment with 
respect to the surface plane. The AHE is absent since we 


consider a non-spinpolarized 2D Rashba electron ga^^^ 
and is not further considered in this paper. Also, we pro¬ 
vide a phenomenological functional form for the different 
components of the residual resistivity as function of the 
orientation of the magnetic moment. 

The paper is organized as follows: In Sec. [H] we briefly 
introduce the Rashba model, basically to define all quan¬ 
tities. Then, the ingoing and outgoing scattering states 
of the Rashba Hamiltonian are introduced. We express 
the resistivity tensor components in terms of scattering 
coefficients, which we relate to the scattering matrix via 
the Lippmann-Schwinger equation. We introduce expres¬ 
sions for cylindrical and non-cylindrical impurity poten¬ 
tials of magnetic impurity moments normal to the sur¬ 
face and arbitrary orientation of the magnetic impurity 
moment. The result of the resistivity tensor i.e. the ap¬ 
plication of the derived expression for an Ee impurity on 
Au(lll) is reported and discussed in Sec. HI. A summary 
is provided in Sec. |IV[ 

II. THEORY 
A. Rashba Model 

The SO interaction leads in a structure-asymmetric en¬ 
vironment such as a surface or interface to a spin-splitting 
of the otherwise two-fold degenerate eigenstates of a two- 
dimensional electron gas. The model of Bychkov and 
Rashbafil^ describes this splitting by adding to the ki¬ 
netic energy of the free electrons the so-called Rashba 
Hamiltonian 

H = (8) I 2 - -j^y^^Py ~ ^yP^) ’ (2) 

where ^ {x, ^}, are the components of the momen¬ 
tum operator p in a cartesian coordinate system with x, y 
coordinates in the surface plane whose surface normal 
points along Cz. m* is the effective mass of the electron, 
cr^ are the Pauli matrices and I 2 is the unit matrix in 
spin-space with a global spin frame of reference where 
the spin z-direction is aligned parallel to Cz. <aso is the 
Rashba parameter, a measure of the strength of the SO 
interaction and the parameter that controls the degree of 
Rashba spin splitting. 

The eigenstates corresponding to this Hamiltonian are 
written as a product of a plane wave in space and a two- 
component spinor 

(3) 

i.e. they can be considered as a superposition of spin- 
up and down-states when measured with respect to 
the surface normal. The orientation of the local 
spin-quantization axis, given by the expectation value 
n±{k) = = ±(-siting, cos= ±e^ 

lies in the surface plane and is perpendicular to the wave 
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vector k = /c(cos 0^, sin 0) = kck. We find that the 
quantization axis is independent of the magnitude k and 
depends only on the direction Ck of the wave vector k. 
With respect to this quantization axis that is parallel to 
60 in a cylindrical coordinate system, are spin pure 
eigenstates and we can associate pJa for = +(—) as 
spin-up (-down) state. 

The energy dispersion is characterized by the /c-linear 
splitting of the free-electron parabolic band dispersion as 
denoted: 


£±{k) = ^ ± ttsofc = ^[{k± k,of - ki] . (4) 

Due to the presence of the ^-inversion broken symme¬ 
try and the SO interaction, the origins of the spin-up 
and -down parabola is shifted by the Rashba or the spin- 
orbit wave vector, respectively, k^o = rrfa^o/li? and the 
dispersion relation describes two bands. 


B. Scattering States 


In order to describe the scattering problem of the 
Rashba electrons at a single impurity, it is convenient to 
exploit the cylindrical symmetry of the Rashba electron 
gas, to introduce the cylindrical coordinates (r, (p) with 
radius r and azimuth (p between the vector r and the x 
axis, and to place the impurity at the origin of the coor¬ 
dinate system. For this purpose it is more convenient to 
express the eigenfunctions of the Rashba Hamiltonian ([^ 
in terms of the cylindrical Bessel functions rather than 
plane waves as 


'0£m± ) 


v/2 


Jm{k±r)e^^* \ ... 


where the state is labeled by the energy 6 , band index 
Of = ±, and the orbital angular quantum number m, 
m G Z, rather than the wave vector k or {k^cp^). The 
wave number for band a is defined by 


H- 

II 

-H 

^S0 5 

y with ku = 

2m* £ 1 u2 

^2 “T AvgQ 


so that kj^ — k- = —2kso holds irrespective of the value of 
e{> 0). We note that 'ipsmai'^) is an eigenvector for the 
2 : component of the total angular momentum operator 
jz = ^ with an eigenvalue jz = and 

is the orbital angular momentum operator. can 

be decomposed into an incident and an outgoing wave: 
'ipsma{r) = ' ^^r cach band 

component, the incoming and outgoing wave functions 
are respectively cylindrical Hankel functions of second 
kind (see Ref.^for a similar derivation). 


= 


2 V 2 


and first kind. 


= 


2 V 2 




(7) 


Their phase factors are chosen such that at large dis¬ 
tances (r ^ 00 ) we can express them as: 


pj 


in 

em± 


(7) 


^JA^Tk±r 


^—ik±r ^imcf) 



(8) 


with and describing 2D cylindrical waves in¬ 

coming toward and outgoing from the origin of the coor¬ 
dinate system, respectively. They are related by 

r4m±G) = Ti ^°G(m+l), 


where T denotes the time reversal operator. 

Now, we introduce a localized impurity for convenience 
placed at the origin of the cylindrical coordinate system 
of the 2D electron gas. We describe the elastic scattering 
of the wave function ^pem± from the impurity with the 
Lippmann-Schwinger equation involving real and spin 
space coordinates: 

+ Jdf'df" Go{f,f',s)t{f',f'')'ipsma{r''), ( 10 ) 

where Go(r, r',6) is the Green function of the Rashba 
electron gas and t(r, r') is the transition matrix (t- 
matrix), related to the impurity potential v(r) via 
the Dyson equation: t(r, r') = v(r) S{f — r') -h 

j df" v(r) Go(r, r",6) t(r",r'). In the asymptotic re¬ 
gion where the impurity potential v(r) vanishes, Eq. 
can be written in a simpler form by using scattering co¬ 
efficients C{ma,m'a'). Then, the wave function of an 
incident electron with quantum state ( 5 , m, a) scatter¬ 
ing elastically from a non-cylindrical impurity potential 
placed at the origin is expressed as 

<Pemair^ = + E \/^ ( 7 )• 

, , V 

m,a' 

( 11 ) 

where the factor ^k^'/ka accounts for the fact that the 
incoming and outgoing waves, and carry 

electron current kM /ka rather than unity due to the rel¬ 
ativistic correction of the velocity operator which will be 
discussed below. Here, the scattering coefficients fulfill 
the unitary condition, 

C{ma,miai) C*{m'a',miai) = 

dmm'^aa' 5 

micvi 




(6) 


( 12 ) 
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Specifically, the diagonal elements of the above equation 
with m = m' and a = a', 

\C{ma,miai)\‘^ = 1, (13) 

micni 

ensure a current conservation. 

For a ^ a' the coefficients give weight to the inter¬ 
band transition during the scattering. For m ^ m' ^ the 
direction of k and thus the total angular momentum com¬ 
ponent of the Rashba electrons changes during the scat¬ 
tering process, and the scattering coefficients refer to the 
amplitude of the intra-band scattering. When v(r) has a 
cylindrical symmetry, i.e. v(r) = v(r), the orbital quan¬ 
tum number m is conserved and C{ma,m'a') simplifies 
to C{ma^m'a')djn,m' ‘ The scattered wave function will 
be a linear combination of the spin-splitted eigenstates 
denoted by the + and — bands. This mixing is due to 
the spin-flip inter-band transitions whose origin is the off- 
diagonal part of the Rashba Green Gq function coming 
from the spin-orbit interaction. 


C. Connection to the scattering matrix 


We present the relation between the scattering coeffi¬ 
cients C{ma,ma') and the t-matrix t(r, r') elements in 
the orbital momentum representation. For this purpose, 
it is convenient to express the Rashba Green function in 
terms of solutions of the Rashba Hamiltonian ( [^ in the 
cylindrical coordinate system presented in Sec. |II B( In 
order to derive the Green function, we fix f' and consider 
Go(r, r',£) to be a function of r. Then, Go(r, r',^) are 
found to be a linear combination of the solutions given 
by Eq. <§ and the out-going solutions (so-called irregu¬ 
lar solutions) given by Eq. 0 . Eurthermore, by taking 
account of the cusp condition of the Green function at 
it is easy to derive. 


r = r 


Go(f,r",e) = 


X 


lE^ 


^ma 

<ma 


i{k^ + /c_) 

;) ^emoL I ^ 






\ /./.out I ^i(m+|)7i 
:/ K'rema \ ^ 


for r > r\ 
for r < r'. 


(14) 


Using Eqs. ( p!4| ) and ( [Tq| ), one arrives at this general ex¬ 
pression for the scattering coefficients: 


C{ma^ m'a') = 


^mm' ^aa' T . 


2kn/> 





k , 


(15) 


D. Residual resistivity tensor 

Using the Kubo linear response formalisnnP^ we can 
show that the components of the resistivity tensor 


measuring the potential drop in direction 7 after applying 
an electric field in direction 7 ' in the DG limit cj ^ 0 are 
given in terms of the scattering solution (fsma by 


P 77 ' 


TVUO 


o 9 9 


-ei-naj){fi-fj) 
i,j (16) 


where 7 , 7 ' G the indices iy stand each for 

(e, m,(a), Uq is the surface electronic density, S denotes 
the area of the surface, e is the electron charge, and 
fi = 0{eF — Si) is the occupation number for the en¬ 
ergy level at T = 0 K. Here, Uq is related to the Eermi 
wave numbers of the two bands, kp^ and kp-^ by 




F+ 




(17) 


and the factor I/S' in Eq. ( [16| ) may be regarded as rep¬ 
resenting the impurity number density, if rii is low 
enough. Also, one needs in Eq. (16) the relativistic ve¬ 
locity operator, z.e., 

( Vx = ^ 

In addition to the prefactor cj, the summation over states 
i and j in Eq. (16) gives rise to another factor cj, since Si 


must satisfy the condition, sp — huj < Si < Sp. In spite of 
this, the right-hand side of Eq. (16) takes a finite limit¬ 
ing value in the limit of cj ^ 0 , since the matrix elements 
{pi \ m'^v^\(pj) evaluated for the asymptotic scattering re¬ 
gion (r ^ oc) diverge as 1/uj in the limit of cj ^ 0 , as 
will be demonstrated in Appendix |A| 

The resistivity tensor is related to the energy dissipa¬ 
tion P in the system per unit time by 


P — ^ ^ P 77 ' Ty Ty 

7,7' 


(18) 


where denotes the 2D current density in the 7 direc¬ 
tion. By inserting the matrix elements of the momentum 
operators in the limit of cj ^ 0 given in Appendix A into 
Eq. (16), one can derive the most general expression for 


the diagonal components of the resistivity tensor. 


P77 ~ 


^TiSnie^ 


^ ^ (^m/,m-|-l A l)^Q 


ma,m' a' 


+ C{ma, la")C{m'a', I + la") 

la" 

C{ma, la")C{m'a\ I — l(a")|l^) 


where kpu = ^m(^f) and the plus and minus signs 
correspond to pxx and pj^ respectively. In deriving 
Eq. (19) we replaced in ( [l^ the sum over discrete 

states by ^ Em=-oo Ea=± ^he scatter- 

ing coefficients C{ma,m'a') are implicitly energy depen¬ 
dent. After integrating over ka these coefficients will be 
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taken at the Fermi level, s-p. In the case of the ab¬ 
sence of the impurity, the resistivity vanishes since the 
scattering coefficients are given in this case simply by 

If we assume that the impurity has no on-site SO cou¬ 
pling and also that it has a magnetic moment perpendic¬ 
ular to the surface, then the t-matrix is diagonal in spin 
space 


t = 


f 0 

V 0 hi 


( 20 ) 


i.e. spin-up and spin-down electrons scatter differently at 
the impurity. However, the scattering is not spin conserv¬ 
ing, because the spin of the Rashba electrons lies in the 
plane and the t-matrix is not diagonal anymore in that 
spin frame of reference. Furthermore, if the impurity 
potential is cylindrical, the orbital momentum represen¬ 
tation of tcra with a G {t, i} reads 


haihr') = 


( 21 ) 


In this case, as seen from Eq. (15), the scattering coef¬ 
ficient C{ma,m'a') becomes diagonal with respect to m 
and m', and the expression for the resistivity, Eq. (19), 
is further simplified as 


P 77 ~ 


S7rSn‘^e‘^ 


E 


a' 

^aa' + C{ma, ma")C {m!a'^ m'a") .(22) 


Obviously, the two diagonal components, pxx and pyy^ 
are identical in this case. 


E. In the limit of feo = 0 

Here, we consider the limit of ^ 0 to derive the ex¬ 
pression of the residual resistivity induced by a localized 
impurity for the 2D free-electron gas without the Rashba 
SO term. This may be useful since the derivation of this 
quantity has not appeared in the literature to our knowl¬ 
edge. Eor this purpose, it is better to choose energy e, 
orbital angular momentum m, and spin index a as the 
quantum numbers for the description of scattering states, 
where the spin quantization axis is chosen as the 2 : axis 
as in previous sections. Namely, instead of Eqs. and 
0, we employ 

(23) 

and 

= ie^-ff^i)(fcr)e*™^|<T), (24) 

as the incident and scattered electron wave functions, 
where k = \/2m*e/h, \ t) = (1, 0), and | 1) = (0,1). 


Then, the wave function of an incident electron with 
quantum state a) scattering elastically from a non- 

cylindrical impurity potential placed at the origin is ex¬ 
pressed as 

f>ema{h) = hZiaih) + E (^ • (25) 

m' ,a' 


By following the same procedure, we can easily show 
that the resistivity tensor for the present case with o^so = 
0 is given exactly in the same form as Eq. ( p!^ , except 
that kFa is replaced by k^ = \/2m*5F/^, the scatter¬ 
ing coefficients of the type C{ma,m'a') are replaced by 
(7(mcr, m'cr'), and further the summation over band in¬ 
dices is replaced by the one over spin indices. 

Eurthermore, if the t-matrix is diagonal with respect to 
electron spin and the impurity potential has cylindrical 
symmetry, we can derive a more simplified expression 
corresponding to Eq. (22), 


P 77 ~ 


hk^ 

^TiSnle^ 


E 


l+(7(mcr, m(j)C''{m!(j^ m'a) 


(26) 


where it should be noted that in contrast to the Rashba 
electrons with a finite ago, the spin-flip scattering does 
not occur in the present case. The scattering coefficient 
in the above equation can be expressed by using the phase 
shift as 


2 


C{ma,ma) = 


(27) 


By inserting this expression into Eq. (26), we obtain fi¬ 
nally 


P 77 


2% 


+00 


E E sin^ [(5„+i(eF,cr) - Sm{eF,cr)], 


cr="|"4, rn= — oo 


(28) 

where we used the relation k^ = 27rne. This is a mod¬ 
ification of Eriedel’s resulfP^ for the residual resistivity 
of a single impurity in a 3D electron gas (see Eq. 0) 
to the case of an impurity in a 2D electron gas with¬ 
out the Rashba-type SO term. The only difference is the 
scattering phase space of momentum transfer in the held 
direction, which is larger in the 3D case than in 2D and 
this is taken care of in Eq. 0 by the multiplicity i 1 
of each angular momentum component. 


F. s-wave scatterer 

In this section we will consider the scattering of Rashba 
electrons by an impurity whose spatial extent is much 
smaller than the Fermi wave length. For such a scatterer, 
one may be able to employ the (5-function approximation 
for the t-matrix, 

h<T'{r,r') = 6{hS{r') 


(29) 
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It should be noted that within this s-wave approximation, 
only '0£,m=o,a(^) having Jo{kar) for its up-spin compo¬ 
nent and 'ips,m=-i,a{'k) having Jo(^a^) for hs down-spin 
component make nonzero contributions to the matrix ele¬ 
ments of the t-matrix, {^srn'a' \ sincc Jm(0) van¬ 

ishes for m > 0. 

We aim at deriving the general expression of the im¬ 
purity resistivity when the t-matrix is given by Eq. (29). 


First, we note that by using Eq. (15), the scattering co¬ 
efficients C{ma^ m'a') for m and nv equal to 0 or —1 are 
given by 


C{0a,0a ) = -daa' --% 


C(y —Icr, —1 q^ ) — H“ 


^/ko, k(y 

2kM 


-s(Q!)s(Q;')iu 


C(0q:,- la') = 


C(-la,0a') = 


(30) 


where s(a) is defined by s(±l) = Tl and the energy 
argument for the t-matrix is omitted for simplicity. For 
m and m' larger than 0 or smaller than —1, we have 

C{ma, m'a') = W5c«a'e-<'"+5)-. (31) 


The expression of the impurity resistivity can be ob¬ 
tained by inserting Eqs. ([so]) and ([H]) into Eq. (p^. 


where twelve pairs of (m, rnj make a non-vanishing con¬ 
tribution to the resistivity. After a lengthy calculation, 
one obtains 


Pjj ~ 


27TSriie^ 




with 


N = Re(f;^f^t) 


-iRe 


(ift - til) ht - (tn - hi) h 


1 2 


tt 


(33) 


(34) 


where the t-matrix should be evaluated at the Fermi en¬ 


ergy and the negative and positive signs in Eq. (32) corre¬ 


spond to pxx and pyy^ respectively. In deriving the above 
equations, we have used the general relation for the t- 
matrix (optical theorem), t — t’l' = —i tt’*', implying in 
the case of a 2x2 matrix that 

Im % = > 

Im ^4,4, = (|f4,4,|2 + , 

th -hf = i ihftn + hiihf) ■ 


Because of the above relations, M and N in Eqs. (33) 


and (34) may be expressed in many apparently different 


but equivalent ways. 


III. RESULTS AND DISCUSSIONS 

As aforementioned, we investigate the example of Fe- 
adatoms on Au(lll) surface with an area S = 1 
(unit area), considering the scattering of the Shockley 
surface states of Au(lll) at the adatom. We assume 
that the impurity has no on-site SO coupling term. Fur¬ 
thermore, we will consider only cylindrical potentials. In 
the case where the magnetic moment of the impurity 
is oriented in the 2: direction, the t-matrix can be ex¬ 
pressed by Eqs. (20) and (IT]). Since the characterisic 
wave lengths of the Rashba electrons at correspond¬ 
ing to and kp- are much larger than the impu¬ 
rity size, we can proceed with the so-called s-wave ap¬ 
proximation (see section HE). Indeed the wavelengths 
Af+ = 27r/0.192 = 32.7 A, Ap- = 27r/0.167 = 37.6 A are 
large considering k^o- = —0.0125 A“^ with parameters 
taken from Reffl The advantages of this approximation 
are the fast numerical evaluation of the scattering coef¬ 
ficients and an easy tracking of the impact of scattering 
on the resistivity. The connection to the phase shift will 
be needed in the up-coming discussion and is given via: 
^crcr(^) = — 1), which are the diagonal elements 

of the t-matrix in spin space (see Eq. 20). 

The s-wave approximation has been used numerous 
times for the interpretation of scanning tunneling mi¬ 
croscopy based experimentThis was done in 
the context of standing waves on Cu(lll) surfac^^ 
or Au(lll) surfac^^ and confined electronic states in 
corrals of Fe or Co adatoms deposited on Cu(lll) 
surfacd^^^. For Fe adatoms on Cu(lll) surface, good 
fits to the experimental features were obtained with a 
phase shift of 7r/2 but a better agreement was found with 
a phase shift of ioo, which would correspond to maximaly 
absorbing adatoms (black dots]P^. In the latter case, the 
overall scattering amplitude reduces by a factor of 2 com¬ 
pared to a phase shift of 7r/2. 

We follow a description similar to the one of Heller et 
al.l^but considering a phase shift of 7r/2 in the minority- 
spin channel. The majority-spin channel is considered 
fully occupied in the magnetic case and thus the cor¬ 
responding phase shift is set to tt. These assumptions 
were used in Ref.l^ and are contorted by our ab-initio 
simulations based on density functional theory as imple¬ 
mented in the Korringa-Kohn-Rostoker Green function 
methocP^. From these calculations we learned that the 
easy axis of the Fe magnetic moment is out-of-plane. The 
adatom local density of states (LDOS) is characterized 
by a resonance close to e-p in the minority-spin channel 
and the exchange splitting between the majority-spin and 
minority-spin resonances is about 2.8 eV. The broaden¬ 
ing of the resonances in the magnetic case is 0.6 eV in 
the majority-spin channel and 0.4 eV in the minority spin 
channel while in the non-magnetic case, it is considered to 
be 0.6 eV. Therefore, we assume that the LDOS consid¬ 
ered in our model follows the Lorentzian shapes depicted 
in FIG.j^a). We note that in our scheme based on linear 
response theory, only the phase shift at the Fermi energy 
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FIG. 1. Local Density of States of an Fe adatom deposited on a Au(lll) surface described by a Lorentzian model wherein 
the broadening is induced by hybridization effects among the electronic states of the impurity with those of the substrate. 
Two cases are considered, a magnetic (a) versus a non magnetic impurity (b). After defining the phase shifts at the Fermi 
energy in the magnetic case, the phase shifts in the non-magnetic case are derived considering the same charge for both type 
of impurities. 


is essential. 


In order to evaluate the impact of magnetism on 
the residual resistivity, we consider a non-magnetic Fe 
adatom and use charge conservation in order to guess the 
appropriate parameters. The spin-dependent charge No- 
of the impurity is given by ^ 6 o{sf^o-) with the Friedel 
sum rule. Charge conservation imposes then that in the 
non-magnetic case = X’ which 

leads to the LDOS plotted in FIG. [^b). 


Considering the approximations mentioned above one 
can investigate the residual resistivity for both cases: 
magnetic and non-magnetic Fe adatoms. To start our 
analysis, we consider a magnetic moment pointing per¬ 
pendicular to the surface. Also to make our study gen¬ 
eral, we explore different SO coupling strengths, which 
then would correspond to the deposition of the impuri¬ 
ties on different substrates. To make this type of inves¬ 
tigations consistent with each other, the energy of the 
highest occupied state of the Rashba electron gas mea¬ 
sured from the bottom of the energy dispersion curve, 
^Rashba _ ^ ig ggt to a Constant, 480 meV 

for the case of Au(lll) surface state characterized by an 
effective mass m* = 0.255mJ^. By changing Sf is 
modified such that ^^^shba vary. 


By inserting —1), — 1), 

and = 0 into Eq. (32), we obtain as the resis¬ 

tivity induced by a magnetic adatom with its magnetic 


moment pointing to the normal direction, 
P77 ~ 




^FM 


sm{6^) sin(^;) cos{6^ - 5i)) 


, (35) 


where 6 ^ and 6 ^ are abbreviations of t) and i 

), respectively. Thus, for the present non-magnetic 
adatom with 


2hk^M 

^7. = ;^ 




(36) 


while for the present magnetic adatom with = ir and 
Si = 7r/2, we simply have 




(37) 


where k^u and Ue are related to ^^^shba^ which is 
kept constant in the numerical calculation, by k^^ = 

and Ue = ^ ^^2 -The latter leads 

to a quadratic decrease of the prefactor ^^^ 2^2 and thus 

of the residual resistivity with respect to k‘^^. This sim¬ 
ply indicates that the more available electrons, the more 
conducting the system becomes. 

The intriguing dependence of Ue on the SO coupling 
strength can be traced back to the particular behavior of 
the density of states of the Rashba electron gas, which 
is characterized by two regimes induced by SO and de¬ 
fined by the two regions of the energy dispersion curve 





















FIG. 2. Evolution of the diagonal components of the resis¬ 
tivity tensor as a function of the spin-orbit wave vector length 
(kso) for a magnetic and a non-magnetic Fe impurity. 


that show a crossing at /c = 0 (see Eq. At en¬ 

ergies below the crossing, the corresponding density of 
states follows a quasi one-dimensional behavior where a 
van Hove singularity occurs at the bottom of the bands. 
Above the crossing, the density of states is a constant as 
expected for a 2D electron gas. By increasing the SO cou¬ 
pling strength, the quasi one-dimensional region becomes 
larger in order to keep constant, which leads to 

the quadratic dependence of Ue on kso and explains the 
strong drop of the residual resistivity when increasing the 
SO coupling strength. 


The latter can be observed in Fig. where the longi¬ 
tudinal residual resistivity as function of SO is depicted. 
The transveral residual resistivity is not shown since it is 
zero for the two cases considered: magnetic (out-of-plane 
moment) and non-magnetic adatoms. Interestingly, mag¬ 
netism and SO coupling strength have opposite impact 
on the residual resistivity. This holds for spin-dependent 
phase shifts that conserve the number of electrons N of 
the impurity after spin-polarization. Indeed, as may be 
seen from Eq. (35), in contrast to magnetism, SO cou¬ 
pling tends to decrease the resistivity. The largest re¬ 
sistivity is found when the SO interaction is switched 
off, which would correspond to the case of Cu(lll) sur¬ 
face. Here the residual resistivity is independent from 
the magnetic nature of the adatom as can be deduced 
from Eq.(35) 


P77 


Sn^e^ 


sin^((5t)+sin^((54,) , 


(38) 


Now, we address the dependence of the residual re¬ 
sistivity on the magnetism of the impurity by analyzing 
the different scattering processes allowed at the Fermi 
surface. The possible elastic scattering processes can be 
found by evaluating the probability for an electron scat¬ 


tering from a state to a state 

/ 27r 

<5(7a - 7'a')’ (39) 

where are given by Eq. (© and a is the band index. 
If the impurity is non-magnetic, the diagonal elements of 
the t-matrix in spin space are equal: = t. In 

this case, the electron scattering probabilities are given 
by 

+ cos (0^-0^,)) 

(40) 

where aa' equals 1 for intra-bad scattering transitions 
(a = a') or -1 for inter-band transitions {a 7^ a'). This 
equation shows that inter-band and intra-band transi¬ 
tions flipping the spin are not allowed since in these cases, 
~ ~ ^ with aa' = 1 and ~ ~ ^ with 

aa' = — 1, respectively (see Fig. 1(a)). 

In the case of a magnetic impurity with a moment per¬ 
pendicular to the surface, the t-matrix is given by Eq. 
and all transitions are allowed, even those leading to a 
spin-flip, as depicted in Fig. (b). Here the electron 
scattering probabilities 

= f l^tt + (41) 

which is different from zero independently from the value 
of (j)^—(j)p . This is due to the magnetic moment of the im¬ 
purity which breaks the time-reversal symmetry. Thus, 
there are more scattering possibilities than in the non¬ 
magnetic case, and therefore for magnetic impurities a 
higher residual resistivity is expected as to non-magnetic 
impurities in-line with Fig. 

Up to now, the magnetic moment was considered 
perpendicular to the surface plane. To generalize our 
study, we explore the impact of arbitrary orientations, 
of the impurity moment, M, on the residual re¬ 
sistivity. The resistivity is a tensor, and contrary to 
the case of an out-of-plane magnetic orientation, its off- 
diagonal elements become finite for arbitrary magneti¬ 
zation directions giving rise to the AMR and PHE. To 
tackle this problem, we rotate the impurity magnetic 
moment pointing initially normal to the surface plane, 
by means of the conventional 3x3 rotation matrices 
IZ G SO(3) by a polar angle 0^ between the direction of 
the magnetic moment and the ^-axis and an azimuthal 
angle (j>^, = 7^^(</>J^)7^J;(6>J^) e^. This translates to 

a unitary transformation of the t-matrix in spin space 
t'(f,r') = U(0^,<^^) t(r,r') Ut(0^,<(,^). U(0^,<^^) = 
U(7^^, (/)j^)U(7^^, are the conventional rotation op¬ 
erators in SU(2) whose representation in terms of a 
2x2 matrix is given by U(7^^,/3) = cos{(3/2) 0 I 2 — 
i sin(/3/2)cr^. The t-matrix for an arbitrary rotation an¬ 
gle of the moment can then be expressed in terms of the 
spin diagonal elements, t^j_and describing the out- 
of-plane moment (see Eq. as: 

t = (8) I 2 + l(ttt - hi) a-e^. (42) 
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FIG. 3. Fermi surfaces scattering processes of Rashba elec¬ 
trons at a non-magnetic impurity (a), magnetic impurity with 
an out-of-plane magnetic moment (b), in-plane magnetic mo¬ 
ment pointing along the a:-direction (c) and along the y- 
direction (d). The transitions between circles with different 
colors are interband transitions, while transitions between cir¬ 
cles with the same color are intraband transitions. The green 
arrows indicate the connection between the initial and final 
state. The cosses indicate prohibited scattering processes, 
while the black arrows at the center of the Fermi contours 
represent the direction of the impurity magnetic moment. 


For given values of matrix t'(r, r') might 

have non-zero off-diagonal components. We note that we 
chose to define the azimuthal angle with respect to the 
x-axis being the direction of the perturbing current. 

The longitudinal and transversal components of the 
residual resistivity in the whole phase space of rota¬ 
tion angles is depicted in FIG. (a) and (b). The im¬ 
ages exhibit a clearly visible angular dependence. The 
anisotropy of the resistivity is in the order of 10“^ nQ. In 
case of the longitudinal resistivity this anisotropy mod¬ 
ifies the isotropic contribution of the longitudinal re¬ 
sistivity, which is in the order of 3.180 x 10“^ nQ by 
about ±0.14 %. Since for the transversal resistivity the 
isotropic contribution is exactly zero, the anisotropy are 
given as absolute values. 

Now we turn to the analysis of the angular dependence 
of the resistivity anisotropy. To the best of our knowl¬ 
edge, no phenomenological functional form for such a 
general angular dependence is available in the literature 
contrary to the case, for example, where the magnetiza¬ 
tion is lying in-plan^^. Therefore to simplify our anal¬ 
ysis, we focus first on the particular orientation of the 
magnetic moment along the x-direction. There, the t- 
matrix given by Eq|g is expressed as 


1 ^ ttt + Hi ~ ^ 

2 V *tt - ^4.4. +^4.4./’ 


(43) 


where and are the upper and lower diagonal 


components of the t-matrix when the magnetic moment 
points along the ^-direction. As we will discuss be¬ 
low, this gives rise to a non-zero off-diagonal contribu¬ 
tion in the resistivity tensor and contributes to the PFIE 
even without spin-orbit contribution at the impurity site. 
When the magnetic moment is in the surface plane {0^ = 
1^), FIG.j^c) and (d) show the behavior of the diagonal 
and off-diagonal components of the resistivity tensor, re¬ 
spectively, as function of the azimuthal angle (j)^. The 
trace of the resistivity tensor is conserved under these 
azimuthal rotations {pxx ± Pyy = constant, V (j)^ while 
the off-diagonal components are related by pxy = pyx- 

The diagonal components of the resistivity tensor 
(FIG. [^c)) can be fitted with the AMR functional 
form given by Thompson et al.l^ p^^ = p±_ ± (py — 
p^) cos^(0j^), where p\\ and p^ define the residual re¬ 
sistivities when the moment is respectively parallel and 
perpendicular to the current. In our particular case, 
the current is pointing along the x-direction, thus p\\ = 

Pxx{(l>M = 0) = Pxx{(l>M = f )• 

{p\\ — P±)^ the maximal value of the AMR, is a pos¬ 
itive quantity as expected for a Rashba electron gaJ^. 
This can be explained by analyzing the different scat¬ 
tering processes on the Fermi surface when the impurity 
magnetic moment is in-plane (see FIGs. [^c) and (d)). 
As done previously, the idea is to evaluate the scattering 
probabilities for an arbitrary rotation of the magnetic 
moment (see Eqj4|). Here we provide the results ob¬ 
tained for scattering processes from = 0 to 0^, = 0 or 
TT. If 0^, — = 0 only inter-band transitions, i.e. a ^ a'^ 

can contribute: 

(44) 

which is zero if the moment points along the p-direction. 
This is the same result obtained for intra-band scatter¬ 
ing probability, when — (j)^ = ii. For the later 

angle configuration, the inter-band scattering probability 
is non-zero independently from the rotation angle of the 
moment: 

(45) 

To summarize, when the magnetic moment points 
along the current direction (x-direction) the back- 
scattering is due to inter-band and intra-band scatter¬ 
ing. However when the magnetic moment is perpendic¬ 
ular to the current direction then the back-scattering is 
only originating from inter-band transitions, which in¬ 
duces a smaller residual resistivity and therefore gives a 
positive maximal value of the AMR, i.e. py > p±. Similar 
scattering processes are possible when the moment points 
along the x-direction or the ^-direction, which explains 
that the resistivities are the same for both magnetic ori¬ 
entations. 

The off-diagonal components of the resistivity tensor 
(FIG. [^d)) could be fitted with the PHE functional 
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FIG. 4. Evolution of the components of the residual resistivity tensor as function of the orientation of the magnetic moment 
in all 47r spatial directions: (a) longitudinal component {pxx)^ (b) transversal component {pxy)- Every point on the sphere 
corresponds to a given orientation of the magnetic moment. Evolution of the residual resistivity tensor components while 
changing only (j)^ when the magnetic moment is pointing in-plane along the x-direction: (c) longitudinal and (d) transversal 
component. Here we plotted the following cases ~ f (black curve) and ~ f curve). 


forirP^ = (p|| — pi_) cos(0j^) sin(0j^) We notice that 
for the considered polar angles {0^ = ~ f)’ PiY 

changes sign when crosses ^ (FIG.IdFd)). This is ac¬ 
companied by a direction switch of the Hall like electric 
field originating from the PHE. The change in the sign 
of p^Y reduces the energy dissipation P given by Eq. [18 
since p^^ is always positive. 

Let us go back to the general case, where the magnetic 
moment points in arbitrary orientations. As mentioned 
earlier, a phenomenological functional form for the resid¬ 
ual resistivity has not been proposed up to now. In Ap¬ 
pendix [^ we derive phenomenological functional forms 
for the residual resistivity tensor and show that the lon¬ 
gitudinal and transversal parts follow a simple angular 
dependence: 

Pxx = P\\ - (PII - P±) sin^(0M) (46) 


Pxy = {p\\ - P±) cos{(f>^) sin(<).j^) sin2(6>j^). (47) 

These equations describe perfectly the angular depen¬ 
dence plotted for instance in FIG. [^c) and (d) with the 


polar angle ~ f (black curve) and = f (red 
curve). 

Alternatively, one may also derive the angular depen¬ 
dence of the diagonal components of the resistivity tensor 
directly from Eq. (32). By applying the aforementioned 
unitary transformations in spin space, Eq. to the t¬ 


matrix and substituting its matrix elements into Eq. (32) 


one yields the diagonal components of the resistivity. In 
the present case, the second terms of M and N, which 
are proportional to the fourth power of tcra', vanish iden¬ 
tically and one obtains 

- sin((5,-) sin((5^) cos((5-|- - (^4,) 


SO \ • 2 

- j — sin 

k^M J 2 


- sin2((5t - 5i) sin^ 0^ sin^ 4>^ 


,(48) 


where similarly to Eq. [^ and are respectively ab¬ 
breviations of (5o(eF,t) Thus, the magni- 
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tude of the AMR, i.e., p\\ — pj_ in Eq. (46) is given by 


hk‘?. 




(49) 


indicating that the maximal value of the AMR occurs 
when the difference between the phase shifts of both spin 
components becomes equal to 7r/2. 


1659/5-1. S.B. acknowledges funding under the DFG- 
SPP 1666 “Topological Insulators: Materials - Funda¬ 
mental Properties - Devices”. The work of H.L was sup¬ 
ported by MEXT KAKENHI No. 25110006 and by JSPS 
KAKENHI No. 24540328. 


IV. SUMMARY 

Using linear response theory, we have derived a formu¬ 
lation of the tensor of the residual electrical resistivity 
for the particular case of a Rashba electron gas scat¬ 
tering at an impurity that can be magnetic and whose 
magnetic moment can point in any arbitrary direction. 
While the obtained form is general, we applied it to the 
case of an Fe impurity deposited at the Au(lll) surface. 
We performed different types of studies and investigated 
the non-trivial impact of the strength of spin-orbit in¬ 
teraction of the substrate, the role of the magnetism of 
the impurity and of the orientation of the magnetic mo¬ 
ment on the diagonal and off-diagonal elements of the 
resistivity tensor. For instance, we found that after scat¬ 
tering, the planar hall effect and an anisotropic magne¬ 
toresistance occur even without incorporating the spin- 
orbit interaction at the impurity site if the orientation 
of the magnetic moment is not perpendicular to the sur¬ 
face. Also an increase of the spin-orbit coupling strength 
induces a dramatic drop of the resistivity, which is re¬ 
lated to a peculiar behavior of the electronic states of the 
Rashba electrons. Magnetism can increase the residual 
resistivity because of the opening of additional scatter¬ 
ing channels, which were prohibited in the non-magnetic 
case. We derive analytically and generalize the usual 
phenomenological functional forms of the angular depen¬ 
dence of the resistivity tensor elements to the cases where 
the magnetization points in arbitrary directions. Finally, 
by switching-off the spin-orbit interaction, we find a sim¬ 
ple formulation of the residual resistivity very close to the 
one given by Friedel for a 3 dimensional electron gaJ^. 

Our numerical results were obtained in the s-wave ap¬ 
proximation involving a Rashba Hamiltonian and as dis¬ 
cussed in the context of lifetime reduction of surface 
states by adatom scatterin^^, it would be interesting to 
investigate the impact of realistic band structures com¬ 
puted from density functional theory on the residual re¬ 
sistivities and asses thereby the effect of other scattering 
channels besides the ones involving only surface states. 
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Appendix A: Evaluation of the momentum operator 
matrix elements 


In this appendix we will calculate the matrix elements 
of the momentum operator between two quantum states 
i = {e^m,a) and j = (e + in the region 

where the scattered wave function can be expressed by 
the asymptotic form Eq. © and extract only the di¬ 
verging terms ^ 1/uj in the limit cj ^ 0. In cylindrical 
coordinates the leading components of the momentum 
operator behaves at large distances as 


h , d h ^ 

Vx -- COS0 ^ H- -kso (Ty 

im^ or m* 

h . , d Ti 

^ ~ T sin 0 7— 

imr or m* 


(Al) 


With this representation of the momentum operators, we 
found that the diverging terms arise from the combina¬ 
tions 


I in,out 

k^s-\-huj,m' ,a' 


m vx Ip, 

^m' ,m 


in,out 


ku 


2iAk 


kn 


(A2) 


where Ak = ka{s huj) — ka{s) ^ 2m*ci;//cM(^) is the 
same for both bands (ct = ±). 


By combining Eq. © and Eq. ( |A2| ) we obtain the 
momentum matrix element in the limit of cj ^ 0, 




Ti kM 

2iAk,/kjpr 


S^{ma, m'a'^ 5), 


(A3) 


where S^{ma^m'a'^e) is given by 

S (jTlOL^m Ct , 5 ) = 4“ 1 )^aa' 

+ C{ma, la")C{m'a', I + la") 

la" 

+ 'y^ C{ma^ la")C*{m!a'^ I — la"). (A4) 

la" 


J.B. thanks M. dos Santos Dias and S.B. thanks Yuriy 
Mokrousov for fruitful discussions. J.B. and S.L. grate¬ 
fully acknowledge funding under HGF YIG Program VH- 
NG-717 (Functional Nanoscale Structure and Probe Sim¬ 
ulation Laboratory-Funsilab) and the DFG project LO 


Analogously, the matrix element for the y component of 
the momentum operator is given by 

(iPj\m*Vy\ipi) ~ (A5) 

2Ak ^ka.k'„ 
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with defined by 

S^irriOL^TIfl , (S) — ^mt ^m — l)^cy.a' 

+ C{ma, la")C*{m'a', I + la") 

la" 

- Y C{ma, la")C*{m'a', I - la"). (A6) 

la" 

Appendix B: Phenomenological derivation of the 
functional forms 

Here we derive phenomenologically the functional 
forms, which fit the computed longitudinal and transver¬ 
sal components of the residual resistivity tensor. The sys¬ 
tem of interest is an adatom with a tilted magnetic mo¬ 
ment interacting with a gas of Rashba electrons. We as¬ 
sume a 2D current density flowing along the x-direction, 
Jxj that generates an electric held E (see Figga)). Be- 
fore analysing the general case of a tilted magnetic mo¬ 
ment, let us recap what is expected when (i) the moment 
lies in-plane and (ii) the moment points out-of-plane. In 
case (i), we proceed as done by Thompson et al.^^ and 
consider the x-component of E: 

cos (p^ + E_i_ sin (Bl) 

where E\\ and E± are the components of the electric held 
parallel and perpendicular to the projection of the unit 
vector of the magnetic moment, lying in the {xy) 
plane (the surface plane) as depicted in FigJ^b). In 
terms of the current density and resistivity, the previ¬ 
ous equation is then rewritten considering the parallel 
and perpendicular projection of the 2D current density 
on the direction of the magnetic moment: 

= J\\P\\ cos(p^ 4- JjLp_i_ sin(/)j^ (B2) 

as function of the azimutal angle ^j^Also, knowing that 
J = J|| cos0j^ = cos0j^ (see Figj^b)) leads to: 

Ei"^ = J {p\\ cos^ <p^ + /9_l sin^ (p^). (B3) 

Here, though, we give this expression in terms of the 
unit vector, defining the orientation of the moment: 

e!^^ = JiPiiie^ • 4)' + P±ieu ■ ey)% (B4) 

Our proposal is that in the general case of a tilted 
magnetic moment, the previous two equations involving 
holds. However, there is a missing contribu¬ 
tion from the out-of-plane component of the magnetic 
moment. In the extreme case (ii), i.e. magnetic moment 
out-of-plane, we have: 

(B5) 

and a simple generalization leads to: 

(B 6 ) 



(b) ^"^ 1- 


FIG. 5. (a) Geometry of the system considered: a magnetic 
moment rotated by a polar angle 6 and azimutal angle 0. 
The current density J is related to the electric field E via 
the resistivity, (b) Decomposition of the electric field and 
the current density parallel and perpendicular to the in-plane 
projection of the unit vector of the magnetic moment 

As deduced from our numerical investigation, py = 
(see FIG.[^(a)). This can be explained from FIG.[^ 
(b) for M II 2 : and FIG. (c) for M || x, where the al¬ 
lowed scattering processes are the same except for the 
interband scattering which flips the spin but does not 
change the direction of k. The latter affects only the 
spin part of the response function, not the residual resis¬ 
tivity (charge part) that we compute. Therefore we get 
the same residual resistivity for M || 2 : and M || x. 

Now we can add up both contributions (i) and (ii) and 
And: 

^ ' ^y) )’ 

which simplifies into: 

Ex = J{p\\ + {pi. -/?||)sin2())j^sin2 6 »j^). (B 8 ) 

A similar approach can be used to derive the functional 
forms for the transversal part of the residual resistivity 
tensor. Here we address the p-component of E and again 
after starting from the form of Thompson et al.l^ for the 
case (i): 

= E\\ sin (p^ - Ej_ cos (p^, (B9) 

we end up with 

H"'’ = (BIO) 

Since there is no transversal resistivity in the case (ii), 
the contribution E^^^ vanishes and we And: 

Ey = J{pa -p||)cos(/)]^sin</>]^sin^ 6 >^). (Bll) 
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